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After a brief introduction to Deep Inelastic Scattering in the Bjorken limit and in the Regge 
Limit we discuss the operator product expansion in terms of non local string operator and in 
terms of Wilson lines. We will show how the high-energy behavior of amplitudes in gauge theories 
can be reformulated in terms of the evolution of Wilson-line operators. In the leading order this 
evolution is governed by the non-linear Balitsky- Kovchegov (BK) equation. In order to see if this 
equation is relevant for existing or future deep inelastic scattering (DIS) accelerators (like Electron 
Ion Collider (EIC) or Large Hadron electron Collider (LHeC)) one needs to know the next-to-leading 
£^ ' order (NLO) corrections. In addition, the NLO corrections define the scale of the running-coupling 

constant in the BK equation and therefore determine the magnitude of the leading-order cross 
sections. In Quantum Chromodynamics (QCD), the next-to-leading order BK equation has both 
' conformal and non-conformal parts. The NLO kernel for the composite operators resolves in a sum 

of the conformal part and the running-coupling part. The QCD and J\f = 4 SYM kernel of the BK 
equation is presented. 
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It is well established that hadrons are made of asymptotically free particles which are quasi on-shell excitations. 
Their confining properties does not allow to measure them directly, and therefore Deep Inelastic Scattering (DIS) 
experiments are performed to study the structure functions of hadrons related to lowest order in QCD to the number 
density constituents. In DIS a lepton with momentum k scatter off a parton inside a hadron with momentum P 

j- through an exchange of a virtual photon with virtuality Q 2 = —q 2 which probes the structure of the target. From the 

■ measurements of the final momentum k! of the scattered lepton one obtains the momentum q = k' — k transferred by 
the virtual photon to the hadron. A suitable frame in which the parton picture is revealed is the infinite-momentum 
frame of the hadron. In this frame we observe different time scale for the interaction of the projectile with the 
constituents of the target and the interaction among the constituents of the target. Furthermore, Lorentz contraction 
t-h . in the longitudinal direction reduces the target into a pancake and time dilation enhances the lifetime of the fluctuating 
partons so that the hadron constituents are effectively frozen during the scattering with the projectile. In the limit in 
. , which the effective mass of the partons is negligible and therefore have small rest mass (partons are asymptotically 
free) the interactions projectile-partons can be considered incoherent. In the Breit frame the longitudinal momentum 
of the virtual photon is zero (g M = (qo, q± 7 0)) and Q 2 = q\. In this case the virtual photon acts as a microscope 
resolving partons in the transverse plane of dimension of the order of -3-. In the limit in which the energy and 

the momentum transfer to the target go to infinity while their ratio is kept fixed (the Bjorken limit), the structure 
functions evolve according to the DGLAP (Dokshitzer, Gribov, Lipatov, Altarelli, Parisi) equation[l| which is an 
evolution equation in Q 2 : the power resolution of the virtual photon probing the target is inversely proportional 
to the momentum transfer; the higher is the momentum transfer the smaller is the size of the parton resolved by 
the virtual photon in the transverse direction. The DGLAP equation, obtained by resumming radiative corrections 
proportional to (a s \nQ 2 ) n for all n > 1, predicts the increase of the parton distributions and describe the evolution 
of the system towards a dilute regime. Besides, although such evolution equation has been successful in explaining 
many features of HERA data like the increase of the parton distribution, in the region of low momentum transfer 
its predictive power is lost and consequently higher twist corrections (1/Q 2 corrections) are not negligible anymore. 
It has been observed in experiments like HI or Zeus at Hera that the structure functions of gluons and sea quarks 
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FIG. 1: Expansion of the T product of two electromagnetic currents in terms of non-local string operators 

grow with decreasing values of the Lorentz invariant variable xb [2| ■ The small- a; b region is obtained in the so called 
Regge limit. When the center of the mass energy square s and the momentum transfer Q 2 are much greater then 
the mass of the hadron, the Bjorken variable is approximately given by xb = jp~^ ~ g ■ The Regge limit is then 

achieved when the parameter s is taken very large while keeping the momentum transfer Q 2 fixed, consequently xb 
decreases. At high energies (Regge limit) contributions proportional to (cv s lns) n for all n > 1 cannot be neglected 
anymore, and therefore they have to be resummed. The equation governing the evolution of the gluon distribution, 
and of the scattering amplitude at high energies is the BFKL|3[ (Balitsky, Fadin, Kuraev, Lipatov) equation. As 
previously mentioned, the DGLAP equation is an evolution equation towards a dilute regime, and it needs higher 
twist corrections for low values of Q 2 . These type of corrections have some theoretical difficulties. On the other hand, 
if we fix the momentum transfer and we let the hadron evolve with rapidity, the parton density increases. Fixing the 
momentum Q 2 means that the virtual photon resolves in the hadron partons with fixed size in the transverse direction 
while their number increases. The BFKL equation provides a nice description of such an increase in density, but it 
fails to describe the evolution of parton density at very high energies: it consider only the emission process, but it does 
not take into account the possible recombination of partons which occurs at such high density. Indeed, the BFKL is a 
linear equation. The recombination phenomena of partons occurring at very high energies are governed by non-linear 
effects (coherency effects) which can be taken into account only by a non-linear equation. Furthermore, in order to 
recover unitarization of the theory, the system eventually evolves towards a saturation region^. Since coherency 
effects at high energies are important, the Breit frame introduced above is not anymore a suitable frame in describing 
DIS processes. We consider the so called dipole frame: we perform a Lorentz boost in a frame in which the hadron 
still carries most of the energy, but the virtual photon has enough energy to split into a quark anti-quark (qq) pair in 
the color singlet state long before the scattering. In this frame the lifetime of the qq fluctuation is much larger then 
the interaction time between this pair and the hadron. The dipole frame is a natural frame work for the description 
of multiple scatterings which as we have seen become important at high energy. We will describe the propagation of 
a qq pair at high energy through a medium generated by color field and we will show that a semi-classical approach 
to the description of the DIS at small-xs leads to a non-linear equation: the Balitsky-Kovchegov (BK) equation[5|,|6|. 

II. HIGH-ENERGY OPERATOR PRODUCT EXPANSION 

Before turning into the discussion of the semi-classical approach to DIS, let us remind ourselves some usual tech- 
niques used in the description of DIS in the Bjorken regime. This will allow us to introduce similar technique for the 
description of DIS in the Regge limit. The Wilson Operator Product Expansion^ (OPE) is a powerful way to study 
the T product of two electromagnetic currents T{j A1 (x)j !y (j/)} at light-like separation x 2 — > 0. The expansion is in 
terms of coefficient functions, pertubatively calculable, and matrix element of operators composed by fundamental 
fields such as quarks and gluons, not perturbatively calculable. An alternative way to the usual OPE is represented 
by the expansion of the T product of two electromagnetic currents in terms of non-local string operators Q. This 
expansion allows one to extract in a more efficient way higher twist corrections to DIS. As in the case of the usual 
Wilson expansion, all the singularities are included in the coefficients in front of such non-local operators; the expan- 
sion in power of (x — yY L of the non local string operator gives back the Wilson OPE in local operators, and in order to 
separate the higher twist effects one has to expand such non-local string operators on the light-cone. Matrix elements 
of such string operators between two hadron state give the structure function of a given stated and their evolution[Toj 
is governed by DGLAP type equation. Furthermore, since this technique is in coordinate space it preserves explicitly 
gauge, Lorentz and conformal invariance. 

As already mentioned above, due to the two different scales of transverse momentum fcj_, in DIS it is natural to 
factorize the amplitude of the process in the product of contributions of hard and soft parts coming from the regions 
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of small and large transverse momenta, respectively. Usually one introduces a factorization scale [i to separate the 
hard and perturbative part from the soft and not perturbative part. The integrals over k\ > /i 2 give the coefficient 
functions in front of light-cone operators while the contributions from k\ < fi 2 give matrix elements of these operators 
normalized at the normalization point /i. The evolution of such matrix elements with respect to the renormalization 
poit /i is the DGLAP evolution equation. The expansion in terms of string operators is diagrammatically given in 
Fig. [n] where dotted lines represent the gauge link 
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Let us now turn our discussion to the case of high-energy (Regge) limit where all the transverse momenta are 
of the same order of magnitude and therefore it is natural to introduce a different factorization scale that is the 
rapidity. Factorization in rapidity means that at high-energy scattering one introduces a rapidity divide rj which 
separate "fast" field from "slow" fields. Thus, the amplitude of the process can be represented as a convolution of 
contributions coming from fields with rapidity 77 < Y (fast field) and contributions coming from fields with rapidity 
7] > Y (slow fields). As in the case of the usual OPE, the integration over the field with rapidity r\ < Y gives us the 
coefficients function while the integrations over the field with rapidity rj > Y are the matrix elements of the operators. 
A general feature of high-energy scattering is that a fast particle moves along its straight-line classical trajectory and 
the only quantum effect is the eikonal phase factor acquired along this propagation path. In QCD, for the fast quark 
or gluon scattering off some target, this eikonal phase factor is a Wilson line - the infinite gauge link ordered along 
the straight line collinear to particle's velocity n^: 



ir>(x±) = Pexpl ig / du n M A»{ un + x± ) > , 
^ J —00 * 



(2) 



Here is the gluon field of the target, x± is the transverse position of the particle which remains unchanged 
throughout the collision, and the index 77 labels the rapidity of the particle. Repeating the above argument for the 
target (moving fast in the spectator's frame) we see that particles with very different rapidity perceive each other as 
Wilson lines and therefore these Wilson-line operators form the convenient effective degrees of freedom in high-energy 
QCD (for a review, see Ref. [HI]). The expansion of the T product of two electromagnetic currents at high-energy 
(Regge limit) is then in terms of Wilson lines 



T{Ux)j u (y)} = Jd 2 Zl d 2 z 2 I™(zi,z 2 )Tr{UlU%} 

r d 2 Zl d 2 z 2 d 2 z 3 I™>(z u z2, Zs)[tr{UZU%}tT{UlU%} - ^fij}] 



We indicate operators by "hat". The diagrammatic version of eq. ([3]) is given in Fig. Ql] the dotted lines represent 
the Wilson lines while the coefficients functions are represented by the propagation of the qq pair in the background 
of shock- wave (represented in the figure by the red strip). 
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FIG. 2: High-energy operator expansion in Wilson lines 
In Ref. [IH such expansion is performed in the M = 4 SYM theory for two BPS-protected currents O 



4 




FIG. 3: Leading-order diagrams for the small- x evolution of color dipole. Wilson lines are denoted by dotted lines. 
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is the composite dipole with the conformal longitudinal cutoff in the next-to-leading order. The appearance of the 
composite operators is due to the loss of conformal invariance of the Wilson line operator in the NLO. Indeed, the 
light-like Wilson lines U(x±) are formally Mobius invariant and consequently the leading-order BK equation is also 
conformal invariant. At NLO the Wilson line operator are divergent and its regularization introduces a dependence 
on the rapidity and conformal symmetry is lost. In order to restore the conformal invariance we redefine the operator 
Tr{[/J C7J! 7 } by adding suitable conterterms. This procedure of finding the dipole with conformally regularized rapidity 
divergence is analogous to the construction of the composite renormalized local operator by adding the appropriate 
counterterms order by order in perturbation theory. 



III. LEADING ORDER EVOLUTION OF COLOR DIPOLES 



Let us consider now the deep inelastic scattering from a hadron at small xb- As previously explained, the virtual 
photon decomposes into a pair of fast quarks moving along straight lines separated by some transverse distance. 
The propagation of this quark- antiquark pair reduces to the "propagator of the color dipole" that is two Wilson 
lines ordered along the direction collinear to quarks' velocity: U(x±)W(y±). The structure function of a hadron is 
proportional to a matrix element of this color dipole operator 

fri(x ± ,y ± ) = 1 - ±-Tr{U«(x ± )U^(y ± )} (5) 

switched between the target's states (N c = 3 for QCD). The gluon parton density is approximately 

x b G{xb^ 2 =Q 2 ) a {p\W{x u Q)\P) 

where r\ = In ■ The energy dependence of the structure function is translated then into the dependence of the color 
dipole on the slope of the Wilson lines determined by the rapidity r\. Thus, the small-x behavior of the structure 



(6) 
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functions is governed by the rapidity evolution of color dipoles [T3|. At relatively high energies and for sufficiently small 
dipoles we can use the leading logarithmic approximation (LLA) where a s <C 1, a s \iixb ~ 1 and get the non- linear 
BK evolution equation for the color dipoles 0,0]: 

±U(x,y) = ^ jd 2 z {x } x z -g_ y)2 [U{x, z) +U{y, z) 

-U{x,y)-U{x,z)U{z,y)} (7) 

The first three terms correspond to the linear BFKL evolution^ and describe the parton emission while the last 
term is responsible for the parton annihilation. For sufficiently high xb the parton emission balances the parton 
annihilation so the partons reach the state of saturation [4| with the characteristic transverse momentum Q s growing 
with energy 1/xb- 



IV. NLO EVOLUTION OF COLOR DIPOLE 

The NLO evolution of color dipole in QCD performed in Ref. [f 4] is not expected to be Mobius invariant due to 
the conformal anomaly leading to dimensional transmutation and running coupling constant. However, the NLO BK 
equation in QCD has an additional term violating Mobius invariance and not related to the conformal anomaly. To 
understand the relation between the high-energy behavior of amplitudes and Mobius invariance of Wilson lines, it is 
instructive to consider the conformally invariant TV" = 4 super Yang-Mils theory. This theory was intensively studied in 
recent years due to the fact that at large coupling constants it is dual to the HB string theory in the AdSs background. 
In the light-cone limit, the contribution of scalar operators to Maldacena- Wilson line[15[ vanishes so one has the usual 
Wilson line constructed from gauge fields and therefore the LLA evolution equation for color dipoles in the M = 4 
SYM has the same form as (|7|). At the NLO level, the contributions from gluino and scalar loops enter the picture. As 
we mentioned above, formally the light-like Wilson lines are Mobius invariant. We also mentioned that the light-like 
Wilson lines are divergent in the longitudinal direction and moreover, it is exactly the evolution equation with respect 
to this longitudinal cutoff which governs the high-energy behavior of amplitudes. At present, it is not known how to 
find the conformally invariant cutoff in the longitudinal direction. When we use the non-invariant cutoff we expect, 
as usual, the invariance to hold in the leading order but to be violated in higher orders in perturbation theory. In the 
calculation of Ref. |12| we restrict the longitudinal momentum of the gluons composing Wilson lines, and with this 
non-invariant cutoff the NLO evolution equation in QCD has extra non-conformal parts not related to the running of 
coupling constant. Similarly, there will be non-conformal parts coming from the longitudinal cutoff of Wilson lines 
in the Af = 4 SYM equation. Here we present the result for the NLO evolution of the color dipole in J\f = 4 SYM ( 
performed in Ref. [12[) in the adjoint representation (we use notations Zij = Zi — zj and {T a )b c = —if abc ) 
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where a is an arbitrary dimensional constant. In fact, a(rj) — ae 71 plays the same role for the rapidity evolution as 
ix 2 for the usual DGLAP evolution: the derivative ^ gives the evolution equation ((SJ). The kernel in the r.h.s. of 

z 2 z 2 z 2 z 2 i— -i 

Eq. (J8J) is obviously Mobius invariant since it depends on two four-point conformal ratios \ 3 24 and \ 2 I 4 . In 12| 
we also demonstrate that Eq. (|SJ) agrees with forward NLO BFKL calculation of Ref. |16|. Let us now present the 
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calculation for the NLO BK kernel in the case of QCD 
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(9) 



where 6 = ^-AT C — |n/. As previously advertised the NLO BK kernel in QCD for the composite operators resolves in 
a sum of the conformal part and the running-coupling part. 



V. CONCLUSIONS 



We have discussed the DIS scattering in the Bjorken Limit and in the Regge limit; we have briefly reviewed the 
standard techniques used to study DIS, that is OPE in local operators and in non-local string operator. This allowed 
us to introduce the OPE at high-energy of the T of two electromagnetic current in terms of Wilson lines operators. We 
have observed the BFKL equation gives us a nice prediction of the increase of the parton density at high energy, but 
it fails at very high energy since it would violate the unitarity condition. The recombination phenomena of partons 
occurring at very high energies are governed by non-linear effects (coherency effects) which can be taken into account 
only by a non-linear equation. Furthermore, we have seen that in order to recover the unitarization of the theory, the 
system eventually evolves towards a saturation region. Then Breit frame introduced above is not anymore a suitable 
frame to describe DIS processes. Instead, we considered the so called dipole frame which is a natural frame work 
for the description of multiple scattering which are relevant at high energy. The evolution of Wilson line operator 
U(x}_ , is the needed non linear equation for the description of the coherency effect appearing at this regime. This 
evolution equation is the BK equation. At the end we have presented the result for the NLO kernel of the BK equation 
in the case of QCD and in the N = 4 SYM theory. In order to restore conformal symmetry we have introduced the 
conformal composite operator. 
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